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SUMMARY 


The  basic  sequential  bearing  only  tracking  problem  for  targets  and 
sonobuoys  by  a moving  monitoring  aircraft  is  examined  here  for  uniqueness 
of  solutions  for  a given  bearing  observation  set.  When  the  monitoring 
aircraft  describes  only  linear  motion,  an  infinite  class  of  ambiguous 
solutions  arises,  the  ambiguity  only  being  resolved  if  prior  knowledge  of 
at  least  one  of  the  parameters  is  obtainable.  However,  if  the  monitoring 


aircraft  purposely  describes  non-linear  motion,  in  general,  except  for 
certain  linearly  restricted  sets  of  target  motion  parameter  values,  which 
may  be  assumed  to  occur  in  practice  with  probability  zero,  no  prior  knowledge 
of  the  target's  parameters  of  motion  is  required  to  obtain  a unique  solution 
for  tracking  from  a given  bearing  observation  set.  If  the  motion  of  the 
monitoring  aircraft  is,  in  particular,  circular,  or  of  a sufficiently  high 
degree  of  non  linearity,  an  unambiguous  solution  is  always  obtained  with 
no  restrictions  whatsoever  on  the  target  motion  parameters. 
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INTRODUCTION 

Sensors  on  the  P-3C  aircraft,  sequentially  in  time,  measure  bearing 
lines  between  the  aircraft  and  a moving  target  emitter.  The  objective  of 
such  measurements  is  to  determine  the  emitter  position  and  velocity,  assum- 
ing for  simplicity,  straight  line  constant  velocity  motion  for  the  target 
or  sonobuoy. 

This  report  discusses  the  type  of  motion  required  of  the  P-3C  aircraft 
in  order  to  track  successfully.  (A  somewhat  analagous  problem  arises  with 
the  use  of  a single  DIFAR  sonobuoy.  The  sonobuoy  is  stationary  and  linear 
motion  is  assumed  for  the  target.) 

CONCLUSIONS 


When  sequential  bearings  are  measured  between  the  P-3C  aircraft  and 
emitter,  a non-linear  (second  degree  or  higher  polynomial  equation  of  motion) 
flight  path,  must  be  flown  by  the  monitoring  aircraft  to  track  the  position 
(and  velocity)  the  the  emitter. 

If  only  a linear  flight  path  is  flown  by  the  P-3C  aircraft,  than  a known 
target  velocity  is  required  for  tracking,  in  order  to  obtain  an  unambiguous 
solution. 

STATEMENT  PROBLEM 

Bearing  only  sensors  emanating  from  a monitoring  (P-3C)  aircraft  are 
used  to  track  a target,  assumed  to  travel  in  straight  line  constant  velocity 
motion.  Given  several  observed  (generally,  in  error)  bearings  between  the 
aircraft  and  target,  is  enough  information  present  to  track? 

ANALYSIS  FOR  LINEAR  MONITORING  AIRCRAFT  MOTION  AND  LINEAR  TARGET  MOTION 

We  show  that,  in  general,  there  is  no  unambiguous  solution  (even  for 
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Che  zero  error  tracking  case).  From  this  it  follows  that  no  correct  regres 
sion  solution  can  exist  for  tracking  with  bearing  errors  present  for  the 
friendly  aircraft.  However,  if  a priori  one  of  the  parameters  of  motion 
of  the  target  is  known,  the  ambiguity  can  be  resolved. 


Consider  first  the  following  figure  illustrating  the  problem  for  the 
zero  error  case: 


FIGURE  1. 

We  assume  for  j =0,  1,  2,...: 

Knouns:  tj  time  of  j Ch  bearing  fix;  0, 

s speed  of  monitoring  aircraft, 

♦ angle  between  x-axis  and  monitoring  aircraft  path, 

9.  angle  between  j successive  bearing  line  and  y-axls 


We  also  let  rx  , r , be  the  x-,  y-  target  positions  at  time  t^ 


It  follows  that 

r r 


(2) 


= v • t . 
Xj  X J 


t 


r = r + v • t. 

yj  yo  y J 


(3) 


x — X • 
tar O ) - i J 

J r v 
Yj  ' yj 


1, ...»  n (n  i 2) 


i 

t 


Using  (1)  and  (2)  in  (3)  yields  immediately 


(4) 


where 


tan^j)  = U 


V + W-t, 
J 


df  a 

U=»  V - 

vdI  r" 

yo 

W1*  v - s-sinjfc/ 

y 


;j  = l,.»n 


I L 
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are  unknown,  since  r , vx , vy  are  unknown. 

yo 

Clearly,  if  (U,  V,  W)  satisfies  the  equations  so  does  (>U,  W,  >W) , 
for  any  real  > 4-  0. 

Thus,  dividing  the  right  hand  side  of  (4)  through  by  U*0,  we 
obtain  with  a little  algebra 


(5)  t cotQ,  = _V_  + _W 

J J u u 

Hence  we  have  the  matrix  relation: 


^ * t.  J i 1, . . , n 


t1  co  tie 


cn  ectOn), 


(n  > 2) 


for  the  zero  error  situation.  Since  trivially  0 = t0  < t^<  . . < t„  jB 
is  n by2  of  full  rank  2,  and  thus  either  (6)  has  no  solution  or  a unique 
solution  for  V/U,  W/U.  (If  X and  X are  two  solutions,  then  BX  = BX 
implies  0 = B • (X-x')  and  premultiplying  by  (BTB  ) • BT  yields  at  once 

x-x'=o.) 


Because  the  zero  error  "observations"  9j,and  hence  tj- cot^jj^j*!,  ..v  n , 
are  generated  by  a physically  realizable  target,  corresponding  to  some  U,V,W, 
and  hence  X,  (6)  will  have  a unique  solution  for  X when  a target  describing 
straight  line  motion  is  present. 


The  solution  X can  be  written 


tx-  cot^j) 
t2’ cot?2) 


since  we  have  perfect  information, 
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(8) 


X = (BtB)_1.  BT2, 


in  a more  redundant  form. 


However,  for  any  such  X,  there  are  infinitely  many  triples  (U,  V,  W) 

that  give  this  ratio,  since  we  have  trivially  U arbitrary,  V = J • V t W 


U • although  ^ and  .H  are  unique. 


Thus,  for  even  the  zero  error  case,  the  solution  set  of  (4)  is: 


(9) 


I j 

I i 

r • 

£•  j 


U arbitrary,  V = U • A,  W = U • B, 

where  X = is  the  unique  solution  of  (§)  given 

in  (7)  or  (8). 


More  generally,  if  0.  is  observed  in  error  as  0;  and  if  we  assume 
Z.  = ttcot  er  Z±  = t±cot  0i,  I - ^ » a ^ j > 

E(llX^  = Z,^na  Cov/fZlX)  = I , a known 

positive  definite  matrix,  not  dependent  on  X,  for  all  X,  then  (4)  is 

replaced  by 


(10) 

and  hence  (5)  becomes 


tan(0j)  = ^ + error^  ; j = n , 

V + W>  tj 


(11) 


t_,«cot(0j)  =»  ^ ^ tj  + error^  ; J = , 


and  analagous  to  (6)  we  have  the  linear  regression  model 


> 


(12) 


o 

2 


B • X + € , where 


— AG  o 

fe  = 2 - 2 = error, 


E ( € ) = 0 , Cov  ( 6 ) = S 1 , 
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1 


Hence,  the  weighted  least  square,  etc.  estimator  of  X is 

a o T -1-1  T _i  o 

(13)  (a)  X (Z)  = (B1  -Z  • B)  1 - B1  Z H 

(b)  E (X  (2  ) | X)  * X, 

(c)  Cov  (X  (§)  | X)  = (BT  • ”1'  B)'1. 

(See  e.g.,  reference  1,  Chapter  4,  for  complete  details.) 


A A A 

Thus,  analogous  to  (9),  the  set  of  'optimal'  estimators  (U,  V,  W)  of 
(U,  V,  W)  in  (10)  is  given  by: 


(14) 


A A A a A A A 

U arbitrary,  V = U • A , W = U • B , where 


A 

X = 


is  given  in  (13a) 


Again,  we  state;  there  is  an  infinite  number  of  triples  (U, 

A A 

to  the  (unique)  ratios  (V/^),  (W/^) . 


V, 


W) 


leading 


Thus,  for  either  the  zero  error  case  in  (4)  or  more  generally  the 
error  case  in  (10),  the  solution  set  of  (U,  V,  W)'s  is  infinite,  and 
unresolvable  ambiguity  for  all  solutions  holds. 


However,  if  one  of  the  parameters  U,  V,  or  W is  known  a priori,  then 
the  ambiguity  is  completely  resolved:  Since  if  V or  W is  known,  a non- 

linear relation  will  always  hold  in  (4)  or  (10)  between  any  function  of 
9j  and  the  unknown  parameters,  we  consider  without  loss  of  generality, 
and  for  simplicity,  the  situation  when  U is  known: 

In  this  case,  (12)  becomes: 
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I 


' 


(15) 

where 

and 


2'  = B • X'  + € ' , 

H'\ 


2 • if 


^ 1 


o . d£  ° , 

U-tj-cottGj) 


X’  JJj  ; E(6  •)  = o,  COV  (6  ') 

x’  = b)”1*  bt  2'”1  a' 


, etc. 


known. 


A 

No  ratios  appear  and  X'  is  unique. 

Thus,  the  optimal  estimator  of  (U,  V,  W)  in  (15)  is 

(16)  (U,  V,  W),  uniquely. 

A similar  unique  solution  for  U,  V,  W holds  for  the  zero  error  case 
when  U is  known. 


ANALYSIS  FDR  NONLINEAR  MONITORING  AIRCRAFT  MOTION  AND  LINEAR  TARGET  MOTION 

We  show  that  in  general  when  the  monitoring  aircraft  purposely  describes 
nonlinear  motion,  the  unknown  target  generating  the  observed  bearing  lines 
is  essentially  uniquely  determined  by  solving  the  bearing  line  equations. 

Any  ambiguity  in  the  general  case  only  occurs  (in  general  with  probability 
zero,  given  all  possible  admissible  values  equally  probable, a priori)  when 
certain  linear  restrictions  hold  for  the  unknown  parameters  of  target 


motion. 

Knnwns : 

• th 

time  of  j bearing  fix; 

(17) 

i± 

°<»  cj  + °(a  cj  + o<3  tj  + a • 0x  (t_,4  ) 

y> 

•if 

£«  0 + b • o„  (t3  ) 
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y-y  positions  of  friendly  aircraft  at  time  t , °^2,  °(  3,  Pi, 

^2,  a,  b all  known.  , a,  or  b may  be  zero,  but  and  (^2  ^ 0. 

if 

Lnnwns:  r - y - initial  (to)  target  position 


Vx  - X - target  velocity  (constant) 


y - target  velocity  (constant) 


Again^r^  = 0,  by  choice  of  x-y  axes.  (See  Figure  2 below.) 


TARGET 

PATH 


a \ Q 


BEARING 
LINE  AT 


(;:)■ 


.Y-AXIS  ' \yi 


ft)  \:i)  f 


MONITORING 
AIRCRAFT  PATH 


FIGURE  2. 

Equations  (2)  and  (3),  with  the  new  interpretation  remain  valid  here. 


We  first  consider  the  zero  error  case. 


Using  (17)  and  (2)  in  (3),  we  obtain  at  once 


9 


NADC-76379-50 


(18) 

tan  , 

where 

u’  = V. 

V*  ^ r. 

W = V 

V'  -i-  W '•  t j ->j* 


J - 1, . ,n 


Clearly,  in  any  physically  realizable  situation  (U',  V',  W')  corresponds 
to  an  actual  target  and  may  take  on  any  set  of  possible  values  (within 
reason)  and  generate  the  "observations'  9j  , j =l,.,n. 


Thus,  we  can  assume  (18)  has  a solution  for  U',  V',  W',  for  the  given 
set  of  0j  's. 

Suppose,  (U' ' , V'',  W'')  is  some  other  solution  to  (18),  and  first 
suppose  U''  ^ U'.  Thus 


(19) 


tan; 


u* '« tj  -xr 
v"tw\<  -yj 


=!,••»  n 


Multiplying  (18)  through  by  its  denominator  and  (19)  by  its  denominator, 
substracting  and  then  resolving  for  tan  (0j),yields 

(20)  tan (0.)  _ (U’  - U”)  • t, 

(V-  - V”)  + (W  - W”).tj 


cj 


<2 

+ fi  ’ t i 

• 

j 

* • » J ^ 

3 

where 

df 

- 

V’  - 

ILL.  fi  = 

V - W" 

But,  the  same  set  of  S^'s  occurs  in  (18),  and  hence, eliminating  tan  9j 
from  (18)  and  (20)  yields 


10 
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OP  + ® * t j 


= ~ ; j 

V'  + W'  • t - y 
•>  J 


Equivalently, 


1 i n • 


(<s+  a • tj  ) • (u;  t.,  -x4  ) = tj  • (V  + w tj  - y4  ) 

for  j - 1,  2, . . , n. 


Collecting  coefficients  of  the  powers  of  the  tj  's  in  the  above 
identify,  using  (17),  yields: 

(23)  (U,-o(1)  • d • tj  + (-«f2-<a  + (U’-c^)  • ®)  • t/ 

+ (-°/3’(5  - o< ^®)  • tj3  - a-0x  (tj4)  -a 

-tj(o<3  tj3  + a-0x  (tj4)  )-S 

= V’  • t;  + (W  - f±)  * t/  - pitj3  - b • tj  0 (t,3) 


for  j = 1,  2, . . . , n 


Equating  coefficients  of  the  powers  of  the  t4’  s in  the  identity  in  (23) 
gives  the  following  overdetermined  system  in  <3  and  B. 


(24)  (a) 


(U'  - o<1).5 


-o(2  -Q  + (U*  -o^)  • fi  - W - 
*(3  +cL  • IB  = ?2 


and, in  addition,  (p  and  B must  satisfy  the  relation 
) a.O  (t/)  •(?+  t j • tj3  + a.O  (tj4)  ) • £ 


b * tj  • 0 (tj  ) , for  J = 1, . . ,n  . 


11 


NADC-76379-50 


We  now  also  suppose  (temporarily)  that  U'^  . Then  solving  in 

(24a)  for  <3  and  then  in  (24b)  for  & and  finally  substituting  into 
(24c)  we  obtain 


(26)  (a) 


0 - — 

u u'  -*j. 

a = H'  - <*?■-  + 

U'  (U'  -o(1) 


3-V>  *2  » CW’  - @ j)  4V> 

U'  -<*  u'  -*  (U’  -O^)2 


U' 


Simplifying  (26)  (c) , we  can  solve  for  V'  : 


V'  = 


<h  • (U>  - °^2  • (w>  ~ tl*)  ' (U’ 

*3  * (U*  -c^)  + *2a 


W',  U'  arbitrary. 


But,  in  general,  (except  with  probability  zero,  given  all  values 
equal)  V'  need  not  be  related  to  U' , W'  and  as  well  o(^,  °^3»  jSj_»  P 2. 

Indeed,  any  possible  triple  of  values  of  (U',  V',  W')  can  occur,  which 
in  turn  with  the  o(j  s and  's  - also  chosen  separately  - generate  the 
' observations'  9j  *S  . 

Equation  (27)  and  its  consequential  linear  restricting  relations  for 
U' , V' , W'  for  non-uniqueness  to  hold  in  (18)  remain  valid  even  when 

*3  = 0. 

If  U'  = * , then  (18)  becomes 


y 1 - °x  (t|2)  • cotfej) 
i : 

yn  - °x  (t*2)  * cot(9nl 


W'  / , 


X 
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an  ordinary  linear  regression  on  V'  and  W' , implying,  analagous  to  (6), 
a unique  solution  for  V'  and  W' . 

Now  suppose  U'  = U' ' . This  implies  that  if  (U'',  V' ' , W'')  is  any 
other  solution, 

(29)  U’-tj  - Xj 

tan  u?-)  = 

Vr  +W'*tj  -y, 

P'-tj  - 

V”  + V tj  - Jf , 


j=l,  2j..,n  , 

and  hence 

(30)  V'  + W'  • tj  = V"  + W"  -tj  ~7j  j 

ior  j = 1,  I,..,  n # 

This  implies  by  equating  coefficients,  immediately  that  V'  = V"  and 
W’  = W' 


Thus,  in  general,  for  all  cases,  (U'',  V'',  W'')  ■ (U',  V',  W'),  and 
equation  (18)  has  a unique  solution. 

Hence,  when  we  purposely  make  the  monitoring  aircraft  motion  nonlinear, 
in  general  a unique  solution  exists  in  equation  (18)  for  that  target 
generating  the  path. 

The  case  where  error  is  present  is  analagous  and  we  can  thus  apply 
nonlinear  regression  or  linearized  regression  techniques  (direct  filtering 
or  Kalman  filtering)  to  the  model. 
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(31) 


U'  • t; 


tan  (Bj)  - 


+ e j 


V'  + W'  • tj 


- y. 

j 


j = 1,...  n 


to  obtain,  in  general,  unique  estimates  of  U',  V',  W'  (at  least  approximately 
optimal) . 


(For  examples  of  these  techniques,  see  references  2 and  3.) 


The  higher  the  degree  in  of  the  monitoring  aircraft  equations  of  motion, 
the  more  restrictive  the  ' sero-probability ' conditions  on  U',  V',  W, 


because  (from  equating  coefficients  of  the  higher  powers  of  t^  ) of  the 


additional  equations 


(32) 


0*3  (1  + °*2  '■&  = ? 


^\n'  + m-l'$  = m-1 


- B 


? m 


to  (24)  (a)  and  (b). 


Thus,  for  example,  if  m = 3 and  hence  3 ^ 0*  it  follows  that 


(33) 

V' 

^ 2 - 2 J 3 

O’  - 

3 ^ 32 

-I’l  + 

^2  ' V'  = & 3 

* 

U’  • *1 

(u'  - y^)2  * 3 
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Equivalently 
(34) 


- f @2  <*2-03  N 

/ 


. (0 -7  ) 


V'  - ( @2. 

"W3 

W’  * £l  - •/  £_2  - C^2_£_3_N)  + @3  ,(U'  -c/1) 

'o(o  O l -O’  ) oi  1 9 


U'  arbitrary 


If  m ">  4,  then  either  (32)  is  a consistent  linear  system  of  equations 
in  <3?  and  (B  or  inconsistent. 


If  the  * s and  s are  chosen  such  that 


(35) 


(*) . 

V 

1 ^3  ^ 

# 

0 

+ ^2 

Uj\ 

4 

\ 

0 / 

l ^ 1 

for  some  real  ^1,  A 2,  then  (32)  is  consistent,  and  thus  (34)  is  the 
extent  to  which  U' , V',  W can  be  restricted  for  nonuniqueness. 


On  the  other  hand,  if  the  ^i's  and  (?i's  are  chosen  such  that  (35) 
does  not  hold,  then  (32)  is  inconsistent  and  no  (3?  and  ® can  satisfy  (32) 
and  hence  we  have  a complete  contradiction  to  the  assumption  of  nonunnique 
ness  of  solutions  of  (18). 


In  particular,  if  the  monitoring  aircraft  flies  a constant  tangential 
velocity  circular  path  with  angular  speed  GJ  beginning  at  the  origin  at 
tQ  ■ 0 ,the  circle  having  a center  (see  Figure  2)  located  on  the  positive 
y-  axis  at  ( ° J , then  eq.  (17)  becomes 


(36) 


Xj  =*  r*sin(CJ-tj) 

y.  * r*(l  - cos(CJ*tj));  J = l,2,..,n 
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Since  the  sine  and  cosine  are  entire  functions  we  have  the  usual 
T-aylor  series  expansions  valid  for  all  t j =l..,n, 

+«*> 


(37) 


E * t j 


>;  - l Pi-tj' 


Asl 


.I  A.  “J 


where 


°^m 


= T' (-l)A+l 

(2A-I)! 


2.A 


°^2A  - O 

Poi  - r-.n/"-* 

fcX;  i 

£u-l  s O 


4 - 1,  2,  3... 


All  the  analysis  in  eq. 's  (17)  - (35)  is  valid  here  with  the  degree 
of  the  path  being  m * + cC  . Thus  (32)  becomes 


(38) 


^•<8 

-(s 

^s-(8 


= fti 

- o 

= 

- o 


But  this  implies 


(39) 


<£ 


. ?2 


a?  * 

o(3 


?4 

°( 5 


■and 


* # # 9 


Ni»,  we  remark,  for  example,  that 
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r » 

-r  • ^3/3‘ 


while 


(41)  ^ 4 = -r  • /u[ 

5 r • Co  5 /s* 

= -5  • OJ 

Hence  we  have  a contradiction  in  (39)  and  thus  (32)  is  an  inconsistent 
system  for  the  circular  motion  here. 

Thus,  no  0?  and  £8  can  satisfy  (32)  and  consequently  eq.  (18)  has  a 
unique  solution  here. 


AJsc  , we  comment  that  if  the  monitoring  aircraft  motion  is  only 
piecewise  described  by  a power  series  in  time,  similar  uniqueness  results 
must  hold  for  the  solution  of  (18) , by  modifying  the  previous  arguments 
beginning  in  eq.  (21). 

Lastly,  we  remark  that  much  care  must  be  exercised  for  the  type  of 
nonlinear  motion  determined  for  Xj  and  y*  in  equation  (31)  for  the 
presence  of  ill-conditioning  in  the  appropriate  approximation  algorithm 


SUMMARY  OF  ANALYSIS 


Target  always  describes  straight  line  constant  velocity  motion. 


1.  LINEAR  MONITORING  AIRCRAFT  MOTION 


If  all  target  parameters  are  unknown  a priori,  then  whether  bearing 


17 


NADC-76379-50 


observation  error  is  present  or  not,  the  target  is  not  uniquely  determinable 
or  uniquely  estimatable  from  the  data.  (See  (9)  for  the  set  of  all 
ambiguous  parameter  values  for  the  zero  error  case  and  (14)  for  the  general 
error  case.)  If  one  of  the  target  parameters  is  known,  a priori,  then  the 
target  motion  is  uniquely  determined.  (For  X- velocity,  V , of  the  target 
known,  see  (15)  and  (16).) 

2.  NONLINEAR  MONITORING  AIRCRAFT  MOTION 

If  all  target  parameters  are  unknown  a priori,  then  for  both  the  zero 
error  and  general  error  cases  for  bearing  observations,  the  target  is 
essentially  uniquely  determinable  (or  estimable) . The  possible  exceptions 
- 'zero  probability  situations,  given  all  things  equally  possible'  - are 
linear  restrictions  on  the  unknown  parameters.  (See  (27)  for  second 
degree  motion  and  the  more  restrictive  (34)  for  third  degree  motion.) 

When  the  monitoring  aircraft  motion  is  of  sufficiently  high  degree 
(m  > 4),  in  general,  there  are  not  even  any  'zero  probability'  linear 
restrictions  on  the  unknown  parameters  which  can  produce  a nonunique 
solution  to  the  bearing  observation  equations  C(18)).  (See  eq.  (35)  and 
following  remarks.) 

For  the  case  of  circular  aircraft  motion  (see  (36)  - (41)),  the  above 
remark  is  valid,  i.e.,  the  target  motion  is  uniquely  determined  from  the 
bearing  observations,  with  no  exceptions. 
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